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ABSTRACT 

We address the formation of 3D nullpoint topologies in the solar corona by combining Hinode/XRT obser- 
vations of a small dynamic limb event, which occurred beside a non-erupting prominence cavity, with a 3D 
zero-/3 MHD simulation. To this end, we model the boundary-driven "kinematic" emergence of a compact, 
intense, and uniformly twisted flux tube into a potential field arcade that overlies a weakly twisted coronal flux 
rope. The expansion of the emerging flux in the corona gives rise to the formation of a nullpoint at the interface 
of the emerging and the pre-existing fields. We unveil a two-step reconnection process at the nullpoint that 
eventually yields the formation of a broad 3D fan-spine configuration above the emerging bipole. The first 
reconnection involves emerging fields and a set of large-scale arcade field lines. It results in the launch of a tor- 
sional MHD wave that propagates along the arcades, and in the formation of a sheared loop system on one side 
of the emerging flux. The second reconnection occurs between these newly formed loops and remote arcade 
fields, and yields the formation of a second loop system on the opposite side of the emerging flux. The two loop 
systems collectively display an anenome pattern that is located below the fan surface. The flux that surrounds 
the inner spine field line of the nullpoint retains a fraction of the emerged twist, while the remaining twist is 
evacuated along the reconnected arcades. The nature and timing of the features which occur in the simulation 
do quantitatively reproduce those observed by XRT in the particular event studied in this paper Moreover, 
the two-step reconnection process suggests a new consistent and generic model for the formation of anemone 
regions in the solar corona. 

Subject headings: Sun: corona - Sun: magnetic fields - MHD - methods: numerical 



1. INTRODUCTION 

Various dynamic phenomena in the solar corona, as for ex- 
ample soft X-ray jets and specific flares, have been associated 
with magnetic reconnection occurring in a three-dimensional 
(3D) magnetic nullpoint topology consisting of a dome-like 
fan separatrix surface located below the nullpoint and a spine 
field line above it (e.g. Lau & Finn 1990; Antiochos 1998; 
[Priest & Forbes 2002). Observational evidence of 3D null- 
point topologies in the corona is prov ided by, e.g., "saddle- 
like" l oop structures (Filippov| |1999| l, "ellipsoidal flare rib- 
bons" ( [Masson et al . 2009), and "anemone " active regions 
within coronal holes (Shibata 'et al.|[T992b| l. The latter are 
characterized by a full or partial ring of radially aligned bright 
loops which connect the opposite polarities of the region and 
the surrounding coronal hole (e.g. Asai et al. 2008 ). Anemone 
active regions are o ften associated with soft X-ray jets (e.g. 
Shibata et al.|1994 1, which are a strong indication of 3D null- 
pomt reconnection occurring in the corona. Further evidence 
for coronal nullpoint topologies comes from subp hotospheric 
sourc e models of the coronal magnetic field (Demoulin et al. 
1993[ l, and from potential and linear f orce-free field extrap- 
olations of flare and jet regions (e.g. [Aulanier et al.||2000 



Fletcher et al. |200 1 [ [Ugarte-Urra et al. |2007[ IMoreno-Insertis 
et al.||2008| l. Consequently, 3D fan-spine configurations are 
mcreasingly used as the initial condition i n numerical simula- 
tions of 3D reconnec tion, jets, and flar es ( Pontin et al.|2007" 



[Pariat et"al. 2009; Masson et al.||2009] l. The reconnection is 
triggered by boundary motions in these simulations. 
The formation of 3D fan-spine configurations in the corona. 



however, has not yet been studied in much detail. In potential 
or near-potential magnetic fields, a 3D nullpoint configuration 
with a fan and a spine naturally occurs if a magnetic bipole is 
embedded into one polarity regio n of a large-scale bipolar am- 
bient field (e.g. |Antiochos|1998] l. Therefore it is expected to 
form when magnetic flux emerges into regions of essentially 
unipolar fields, such as a coronal hole. The fan-spine con- 
figuration then results from the relaxation of the coronal field 
after its reconnection with the emerging flux. 2D numerical 
simulations of coron al soft X-ray jets , inspired by the flux 
emergence model by Heyvaerts et al. (1977), have demon- 
strated that reconnection between emergmg bipolar flux and 
a vertical or oblique coronal field yields the formation of hot 
loops connecting the ambient fie ld with the opposite polarity 
flux of the e merging bipole (e.g. Yokoyama & Shibata|1996 
INish izuka et al. 2008 ). 

Moreno-Insertis et al. (2008 ) recently performed a 3D sim- 
ulation of the emergence of a twisted flux tube into an oblique 
unipolar coronal field. As in the 2D cases, they found the 
launch of a jet and the formation of a growing system of 
hot reconnected loops connecting the ambient field with the 
emerging flux of opposite polarity. The resulting fan surface 
extends on one side of the emerging region, while on the other 
side it consists of non-reconnected emerged loops only. The 
latter are not strongly heated and would hence unlikely be 
seen in soft X-ray observations, therefore only one half of an 
a nemone loo p pattern should be visible. 



Pariat et al. (2009 ) used an alternative approach to produce 
a jet, by starting from a 3D nullpoint topology and driving the 
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jet by reconnection between open and closed field lines, after 
the latter have been significantly twisted by line-tied boundary 
motions. One outcome of this calculation is that, as a result 
of reconnection of twisted fields, the nullpoint moved around 
the axis of the spine, thus allowing reconnection of field lines 
from all sides of the fan. While this evolution may allow for 
the brightening of the whole fan in soft X-rays, it still does 
not explain how the fan-spine topology was formed in the first 
place. 

These models, although nicely reproducing coronal jets, 
their associated inverse-Y "Eiffel tower" shape, and the field 
line geometry obtaine d from a linear force-free fie ld extrapo- 
lation of a jet region ( Moreno-Insertis et al.|2008] l, do there- 
fore not yet provide a satisfying scenario for the formation of 
anemone active regions. As mentioned above, these regions 
typically show a ring of bright loops below the jet that is remi- 
niscent of a fan dome, which in some cases might well extend 
the area of emerged flux. 

A two-dimensional effect may also play a role in the bright- 
ening of both sides of fan surfaces which form during flux 
emergence. This effect is the reconnection of emerging field 
lines back and forth with the ambient fields on both sides 
of the emerging flux. This is exactly what happens in the 
2.5D simulation of Murray et al. (2009 ;), where an "osc illa- 
tory reconnection" pattern (C raig & McClymont||1991 



curs, which the authors attributed to thermal pressure effects 
around the reconnection layer This process was also found in 
2.5D simulations of a quadrupolar closed field configuration, 
being driven by a non line-tied chromospheric ad-hoc mono- 
tonic force ( Karpen et al.|1998] ). Reconnection back and forth 
was there attributed to coronal relaxation, as a response to an 
"overshoot" due to a chromospheric driving which was faster 
in the simulat i on tha n on the real Sun. In the simulation of 
[Murray et al.| (|2009), the photospheric emergence velocities 
were small compared to the coronal Alfven speeds (Murray, 
private communication, 2009), so that such an "overshoot" 
probably did not occur However, the existence of (non-wave 
driven) oscillatory coronal reconnection in fully 3D configu- 
rations is yet to be established. 

An interesting observational f eature of coronal jets is their 
frequent transverse oscillation (Cirtain et al. 2007 1. Even 
though 2D os cillatory reconnecti on could account for such 
perturbations ( Murray et al.[2009) l, they could also be caused 
by non-steady reconnection m a turbulent current sheet, where 
magnetic islands are gradually formed and destroyed (e.g. 



Yokoyama & Shibata|1996HKhem et al.|2 000; Arc hontis et al. 
2006) , as well as by upward propagating Alfve n waves be- 
ing launched from the reconnection point ( Cirtain et al.|2007 



Moreno-Insertis et al.||2008 1, or by the propagation of a tor 



sional Alfven wave resulting from the reconnecti on of kink 
ing twisted field lines with their ambient field (Pariat et al. 
^09, Filippov e t al. || 20() 9]l. The latter mechanism provides 
an explanation also for the frequently observed helical pat- 



terns traveling along jets (Shibata et al.„1992a{ Patsourakos 
let al.|2008[|Nistico et al.|20n^ 

In this paper, we propose a fully 3D two-step reconnection 
model for the formation of broad fan-spine configurations in 
the corona. The model results from a zero-/? line-tied MHD 
simulation, in which the evolution of the coronal magnetic 
field is driven by twisted flux emergence prescribed at the 
photospheric boundary. The simulation was initially devel- 
oped for the interpretation of a puzzling event observed at 
the solar limb by the Hinode X-ray telescope (XRT). In this 
event, two distinct small coronal loop systems developed one 




Fig. 1 . — (a): MSDP Ha observation of the prominence on 2007 April 24, 
at 13:38 UT. (b): XRT observation of the bright loop systems on 2007 April 
24, at 18:35 UT (compare with Fig.[2l, overlaid with Ha intensity contours 
of (a). 

after the other beside the edge of a prominence cavity, the first 
one apparently "feeding" the second one, while a swaying jet- 
like brightening was propagating along the cavity edge. Our 
model does not only reproduce the shape and timing of the 
main features observed in this particular event, but also ac- 
counts for the formation of full (and not half) anemone bright 
features, as generally observed in soft X-rays within coro- 
nal holes. It finds (in accordance with some past observa- 
tions and models of coronal jets) that a large fraction of the 
emerged magnetic twist reconnects and is evacuated upward 
in the form of a torsional Alfven wave. It furthermore shows 
that nullpoint reco nnection can be accompanied by slipping 
reconnection (e.g. lA ulanier et al.| [2006), which is supported 
by apparently slipping cavity loops observed by XRT in the 
event. It finally predicts that a fraction of the twist eventually 
remains around the inner spine beneath the fan/anemone sur- 
face, which therefore does not fully relax to a potential field 
configuration, even though it looks potential at first order. 

2. OBSERVATIONS 

2.1. Instrumentation 

The observations presented here were obtained during a 
coordinated campaign of prominence observations, involving 
several space and ground-based instruments. They were per- 
formed in the JOP 178 frame on 2007 April 23-29, during the 
first SAJMEK-Hinode observing campaign. JOP 178 has been 



running successfully many times in the past (see http : // 
bass2000 .bagn . obs-mip .fr/jopl78/). A promi- 
nence surrounded by a cavity on the west limb, at S30 
degree, was extensively studied during the campaign (e.g. 
Heinzel et a l. 2008). Fig. [T] shows an Hq observation of the 
prominence obtained by the Multi channel subtractive dou- 
ble pass spectrograph (MSDP) operating at the Solar Tower 
of Meudon, as well as an overlay of the corresponding in- 
tensity contours with the XRT observations described below. 
The prominence was located in a quiet Sun region, apparently 
along the polarity inversion line (PIL) of an extended bipo- 
lar area of weak magnetic field, without visible strong field 
concentrations. A corresponding filament channel could be 
observed days before in EUV. Here we focus on XRT obser- 
vations of a small dynamic event which took place at the edge 
of the cavity on 24 April. 

T he Hinode mission h as been operating since 2006 Octo- 
ber (Kosugi et al. 2007| l. XRT is a high resolution grazing 
incidence telescope which consists of X-ray and visible light 
optics and a 2k x 2k CCD camera. A set of filters and a 
broad range of exposure times enable the telescope to ob- 
serve hot plasma in the range 10^-10® K (for more details 
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Fig. 2. — Snapshots from the XRT observations (300" X 212") at different times. The bright loop systems are denoted by LI, L2, and K. Asterisks and 
horizontal arrows are plotted at the same positions in all images to outline the oscillatory transverse motion of K and the displacement of its elbow. The vertical 
arrow points at loops which are seen to "slip" toward the prominence at the later phase of the evolution. See the text for details. (An mpeg animation of this 
figure is available in the online journal.) The animation uses a different intensity scaling, which outlines the cavity loops in more claiity. 



see Golub et al.|2007 1. The observations presented here were 
obtained with a resolution of 512x512 pixels, each pixel hav- 
ing a size of 1.03" x 1.03". The filter combination used was 
ALpoly/Open. The exposure time was 8.19 s or 16.38 s, the 
cadence was 60 s. 

2.2. Event description 

The brightening of several loop systems at the edge of the 
prominence cavity was observed by XRT for about two hours 
on 2007 April 24, between about 18:30 and 20:30 UT. The 
dynamic evolution described in the following lasted about 45 
minutes, ending at 19:15 UT. 

Fig. |2] shows snapshots outlining the main features of the 
evolution. At about 18:30 UT, two bright loop systems be- 
come visible, one set of small closed loops at the left cav- 
ity edge (denoted as LI), and one large loop arching over 
the prominence (denoted as K), outlining the cavity edge 
(Fig. l2k). K exhibits an elbow (indicated by the horizontal 
arrowjC apparently located above LI. A very small loop sys- 
tem is observed at the apparent footpoint of K (denoted as L2 
in Fig. [2^). In the following minutes, LI continuously grows 
and the elbow slowly moves to the left. At about 18:39 UT, a 
brightening begins to propagate along K, starting from below 
the elbow. 

At about 18:45 UT, the evolution becomes more dynamic: 
LI seems to expand toward L2, loops apparently connecting 
LI and L2 become visible, and L2 starts to grow (Fig. [2]3,c). 
At the same time, the propagating brightening has arrived at 
the elbow and now moves further along K, exhibiting a jet- 
like appearance. The upper part of K moves leftward (see 
asterisks) and after a while the elbow cannot be observed any- 
more. Shortly after, the part of K on the opposite side of the 
prominence starts to brighten weakly, indicating that the loop 
has been filled with hot plasma, ejected from below the elbow 
(the propagation of the jet-like brightening along K is better 
visible in the online movie accompanying Fig.[2]i. 

At 18:52 UT, LI and L2 are roughly of the same size. They 



appear to have a common footpoint region and they collec- 
tively exhibit an anemone-like shape (Fig. |2}l,e). About this 
time, the upper part of K moves back toward the right. At 
19:01 UT LI has started to fade away, and L2 has stopped 
growing (Fig.|2^). From about 19:10 UT on, LI is no longer 
visible. At 19^9 UT, the upper part of K has returned to its 
initial position (Fig.|2p. The transverse oscillatory m otion of 
K is sug gestive of an Alfven wave traveling along it ( Cirtain 
let al.|2007| ). Between about 19:00 UT and 19:20 UT, a suc- 
cessive leftward displacement of loops starting from the foot- 
point of K on the opposite edge of the cavity is observed (see 
the vertical arrow in Fig.l2p. After that, L2 slowly fades away 
in about an hour and no further significant dynamic evolution 
is observed. 

The observations suggest that the dynamic evolution was 
caused by the interaction of newly emerging flux with the 
arcade-like field overlying the prominence. The inspection 
of SoHO/MDI magnetogram data from the days leading up 
to the event did not reveal significant long-lived bipolar field 
concentrations in the vicinity of the prominence, which would 
have indicated the presence of a 3D nullpoint topology before 
the event (as in the simulation of Pariat et al. 2009). Since 
the prominence-cavity system was located at the limb during 
our event, direct observations of emerging flux are not avail- 
able. However, the brightening of the elbowed loop system 
K, as well as the jet-like brightening propagating along it, can 
be understood by means of models of magnetic reconnection 
between a small emerging bipole and a predominantly ver- 
tical coronal field (see Sect. [T]). This suggests that the loop 
system LI was not outlining emerging fields, but rather a re- 
connected arcade that formed along with the elbow in K. The 
fact that LI has a significant height when it becomes visible 
supports this interpretation. Still, the transverse oscillation of 
K, the growth of the loop system L2, and the "slipping"-like 
motion of loops at the opposite edge of the cavity cannot be 
understood straightforwardly within this scenario. 

We note that the observed dynamics did not seem to have a 
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noticeable effect on the stability of the prominence-ca vity sys 
tern. The pr ominence was still observed on April 25 ( Heinzel 
et al.||2008| and on April 26 by the Meudon Solar Tower It 



hence appears that this is a case where emerging flux in the 
vicinity of a filament or prominence does not result in the 
eruption of the latter (see, e.g., ,Feynman & Martin)1995j ). 

3. NUMERICAL SIMULATION 

In order to understand the full dynamics observed, we per- 
form a 3D MHD simulation of the interaction of emerging 
flux with an arcade-like potential field overlying a coronal flux 
rope. The choice of such a coronal topology is supported by 
the presence of the prominence-cavity system. Magnetic flux 
ropes have b een successfully used to model prominences and 
cavities (e.g. Low & Hundhausen 1 995 ; ' Aulanier & Demoulin 
p98 , van Ba ll egooij en 2004; Gibson & Fan 2006|l. We use 
the analytical flux rope model by |Titov & D emoulin ( 1999 
hereafter TD) as the initial condition in the simulation. ITie 
model consists of a toroidal current ring with major radius 
R and minor radius a, which is partly submerged below a 
"photosphere" and is held in equilibrium by an external po- 
tential field created by two subphotospheric magnetic charges 
±q, which are placed at the axis of the torus, at distances ±L 
from the current ring. The coronal part of the model yields 
an arched, line-tied, and twisted flux rope which is embed- 
ded into an arcade-like potential field (see Fig. 2 in TD). The 
depth d of the torus axis (and hence of the magnetic charges) 
below the photospheric layer determines the compactness and 
strength of the magnetic flux distribution in the photospheric 
plane. Here we choose a relatively large depth (see below), in 
order to account for the observed extended area of weak field 
above which the promin ence was lo cated (see Se ct.|2.1 

Previous simulations (|To rok et al. 2004, Toro k &" 

|2005^ '2007';Schrijver et al. 2008) and analytical calculations 
Psenberg & Forbes 2007) have shown that the TD flux rope 
can be subject to the ideal MHD helical kink and torus insta- 
bilities. Therefore we use a weakly twisted rope here, with the 
field lines winding on average only once about the rope axis 
(in a right-handed sense), and we choose the potential field 
such that the rop e is stabl e with respect to the torus instability 
(jBateman 1978; Kliem& Torok 2006). 

As in these previous simulations we integrate the P — 
compressible ideal MHD equations: 



dtp^~V ■ (pu) , 
pdtu — —p ( M • V ) M - 



jxB + W -T, 



(1) 
(2) 
(3) 



where B, u, and p are the magnetic field, velocity, and mass 
density, respectively. The current density is given by j = 
p,Q^ VxB. The term T is the viscous stress tensor, included 
to improve the numerical stability (Torok & Kliem 2003). We 
neglect thermal pressure and gravity since we are interested in 
the qualitative evolution of the magnetic field only. 

The MHD equations are normalized by quantities derived 
from a characteristic length, taken here to be the initial apex 
height of the TD flux rope axis above the photospheric plane, 
R — d, the initial magnetic field strength, Bq, and Alfven ve- 
locity. Van, at the apex of the axis, and derived quantities. We 
use a nonuniform cartesian grid of size [—4,4] x [—5,5] x 
[0,8], resolved by 261 x 301 x 208 grid points (including two 
layers of ghost points in each direction which are used to im- 
plement the boundary conditions), with a resolution of Ax = 
Ay = 2Az = 0.02 in the box center. The resolution is nearly 



constant in the subdomain [—1, 1] x [—1.5, 1.5] x [0, 1], and 
increases exponentially toward the boundaries, to maximum 
values Axmax = 0.14, Ay,„ax = 0.10, and Az^ax = 0.40. 
The plane {z = 0} corresponds to the photosphere. 

The TD flux rope axis is oriented along the y direction, with 
the positive polarity rope footpoint rooted in the half -plane 
{y > 0} (see Fig. pll. The normalized geometrical flux rope 
parameters used are: R — 2.75, a ~ 1., L = 2.5, and d = 
1.75. The top left panel in Fig. [3] essentially shows the initial 
TD configuration, except for the small parasitic bipole and the 
blue field lines on the left-hand side of the TD flux rope. 

We employ a modified two-step Lax-Wendroff method for 
the integration, and we additionally stabilize the calculation 
by artificial smoo t hing of all integration variables (see Sato 
& Hayashi||T979l [Torok & Kliem|[2003| for details). TEe 



reconnection occurrmg m our simulation is due to the re- 
sulting numerical diffusion. The initial density distribution 
is Pq{z) — 2.6exp(— [z + Az]/l.l), chosen such that the 
Alfven velocity, Va, slowly decreases with height above the 
TD flux rope. The system is at rest at i = 0. 

We first perform a numerical relaxation of the system for 
37 Alfven times and reset the time to zero afterwards. We 
then model the emergence of a second twisted flux rope in 
the vicinit y of the TP rope, fo llowing the boundary-driven 
method by Fan & Gibson ( 2003 hereafter FG). In their model, 
a toroidal twisted flux rope is ngidly emerged from a fictitious 
solar interior into a coronal magnetic field by successively 
changing the boundary conditions in the photospheric layer 
of the simulation box. We refer to Fig. 1 in FG for a sketch 
of the model. In our simulation, we choose the FG torus to 
be about one order of magnitude smaller in size than the TD 
torus (the major and minor radius of the FG torus are 0.3 and 
0.2, respectively), in order to account (within the limitations 
given by the finite number of grid points) for the typically 
large difference in size between quiescent filaments and small 
bipole s that emerge in their vicinity (see, e.g., Feynman & 
Martin|I995), which was also suggested by the relative scale 



lengths oTlh e coronal loops observed in the event described 
in Sect.lZ2l 

The FG flux rope is uniformly twisted along its cross- 
section. The twist is chosen right-handed, with the field 
lines winding ~ 4.5 times along the whole torus. We posi- 
tion the rope such that the emergence region is centered at 
{x, y) = (—1.0, —0.5), within the large-scale negative polar- 
ity of the TD potential field. As the TD flux rope, it is ori- 
ented along the y direction, but with the opposite orientation 
of the axial magnetic field (see Fig. |3]l. The magnetic field 
strength within the FG torus varies along its cross-section, 
being ^ 3 So at the outer torus surface, 4Bq at the axis, and 
^ 13i?o at the inner torus surface, for the parameters used 
in the simulation. Although the emergence is driven until the 
apex of the inner surface approaches the bottom boundary of 
the box (see below), the field strength which effectively enters 
the corona during the simulation does typically not become 
larger than ~ GBq. The field strength of the large-scale TD 
field in the small volume above the emergence area is approx- 
imately constant, ~ 0.7 Bq. Within the TD flux rope, the field 
strengths vary between ^ 0.6 and 1.0 Bo- We discuss the ra- 
tionale for our choice of the orientation and strength of the 
magnetic field within the FG flux rope in Sect. [5] 

The boundary-driven emergence is imposeain the layer 
{z = — Az}. Within the emergence area in this layer, we 
overwrite the pre-existing TD field by the respective FG flux 
rope field, and we set the vertical velocity equal to the respec- 
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tive driving velocity, while keeping the horizontal velocities 
at zero. Outside this area, the TD field and the density in 
{z = —Az} are kept at their initial values, and the velocities 
are set to zero, at all times. These settings lead to significant 
jumps in strength and orientation of the magnetic field (i.e. to 
the formation of large values of V • B) at the interface be- 
tween the TD and FG fields at and close to the bottom plane. 
Since our code does not conserve V • B = to rounding er- 



ror, we use a diffusive V • B cleaner ( Keppens et al.|2003 , as 
well as Powell's source term method (Gombo si et al.| 1994 1, to 
minimize unphysical effects resulting from V • B errors. Fur- 
thermore, our overspecified boundary conditions (see above) 
trigger spurious oscillations, which after some time lead to 
numerical instabilities close to the bottom plane, in particu- 
lar at the interface between the TD and FG fields. In order to 
prevent these instabilities, we apply an enh anced smoothing 
of all variables close to the boundary (as in Torok & Kliem 
and we set the Lorentz force densities at {z = 0| to 



2003 



zero at all times. We find that these settings result in the for- 
mation of a rising twisted flux tube above the emergence area, 
as desired. 

The emergence is driven quasi-statically with a maximum 
velocity of O.OlWao- The driving velocity is linearly increased 
and decreased for 10 Alfven times before and after the main 
emergence phase (which lasts for 30 Alfven times), respec- 
tively. The emergence is stopped at < = 50, when the apex 
of the inner surface of the FG torus has reached the layer 
{z = — Az}. The total twist transported into the corona by 
the emerged section of the FG torus corresponds to ^ 1.7 
field line turns. 

4. SIMULATION RESULTS AND COMPARISON WITH XRT 
OBSERVATIONS 

In this section, we describe our simulation results and com- 



pare them with the XRT observations described in Sect. 2.2 
The interaction between the emerging FG flux rope and the 
pre-existing TD coronal field results in t wo di stin ct re connec- 
tion phases, wh ich a re described in Sects. 4.1 and |4.2| respec- 
tively. In Sect. |4.3[ we discuss the magneti c fiel d geometry 
resulting from the reconnections, and in Sect. |4.4| we describe 
the response of the TD flux rope and its surrounding arcade 
to the dynamics accompanying the reconnection. In order to 
emphasize that the evolution found in the simulation does not 
rely on the specific flux rope models used, we refer to the FG 
rope as "emerging flux rope", and to the TD rope as "promi- 
nence flux rope" throughout this section. 

Figs. [3j|5] display magnetic field lines which outline the 
main features of the magnetic field evolution in the coronal 
domain. Figs.[3]and|4]focus on the dynamics, showing the first 
and second reconnection phase, respectively. Fig.lSlshows the 
evolution of the magnetic topology. The blue field lines in 
Figs.[3]-^are integrated starting from the negativepolarity of 
the emerging flux rope, green field lines in Figs. [3]and|4] out- 
line the core of the prominence flux rope (black circles mark 
field line dips). Red (pink) field lines in Figs.[3]and|4]show the 
outer (inner) arcade overlying the prominence flux rope ini- 
tially, and reconnected field lines later on. Note that all field 
lines are calculated from the same positions on the left-hand 
side of the prominence flux rope in all panels. 

4.1. First step: formation of one half of the anemone and jet 
acceleration 

As the emerging flux rope (closed blue field lines in the 
bottom panels of Fig. |3]l enters the coronal domain, it starts 



to expand and a current sheet forms above the outer edge of 
its positive polarity, at the interface of the rope and the outer 
coronal arcade that surrounds the prominence flux rope. Since 
the outermost emerging fields and the outer arcade fields are 
oppositely directed at the location of the current sheet, the two 
flux systems readily start to reconnect, forming a new small 
loop system below the current sheet and strongly bent, elbow- 
shaped field lines above it (cusp-shaped red field lines and 
blue field lines in the central panels of Fig. |3] respectively). 
Note that the not yet reconnected part of the emerging flux 
rope continues to expand in the corona after the reconnection 
has started. The shape of the reconnected field lines agrees 
very well with the shape of the bright loops observed by XRT 
in the early phase of the event described in Sect. 2.2 (see LI 
and K in Fig. [2^), which indicates that the observed loops have 
been formed in an analogous reconnection process. 

The field line shapes are a signature of a fan-spine configu- 
ration in a 3D nullpoint topology (see also Fig. S]?). A mag- 
netic nullpoint is indeed formed within the current sheet in the 
simulation, right after the reconnection has started. It forms 
as the system tends to relax to a lower energy state'. The 
reconnection continues as time evolves, thus the size of the 
reconnected red loop system increases (Fig.|3]l, in agreement 
with the observed growth of LI (Fig.[2^,b). The ongoing ex- 
pansion of the emerging flux rope initially slowly pushes the 
nullpoint away from the prominence flux rope, in agreement 
with the observed leftward displacement of the elbow indi- 
cated by the horizontal arrow in Fig. [2] 

As the emergence continues, the emerging flux rope field 
lines become increasingly sheared with respect to the sur- 
rounding coronal arcade. This is due to the fact that the twist 
within the FG flux rope is nearly uniform along its cross- 
section (as in the well-known Gold-Hoyle model). While field 
lines far away from the flux rope axis are strongly inclined 
with respect to the axis, field lines close to the axis are almost 
aligned with it. Hence, as the flux rope emerges, its outer 
field lines resemble a nearly potential coronal arcade that is 
oriented almost orthogonal to the local PIL, whereas its in- 
ner field lines (i.e. those close to its axis) resemble a small 
sheared coronal arcade (see Fig. [6]). The first flux rope field 
lines which reconnect with the large-scale coronal arcade are 
thus almost unsheared with respect to the arcade (see the bot- 
tom left panel in Fig. [3]l. As the evolution continues, pro- 
gressively more sheared loops are reconnected. As a result, 
the system of new reconnected (red) field lines eventually de- 
velops a shear distribution that is opposite to the one of the 
emerging flux rope: the field lines are sheared at the edges of 
the system, and almost unsheared close to the local PIL (see 
the bottom right panel in Fig. [3] and the corresponding online 
animation). 

The reconnection does not only yield the transfer of twist 
from the emerging flux rope into the newly formed red loop 
system. Part of the flux rope twist is also transferred into the 
lower parts of the reconnected overlying blue field lines that 
are now rooted in the negative polarity of the rope. Since the 
upper parts of these field lines are nearly potential, whereas 
their lower parts experience a sudden injection of twist, they 
are far from being force-free. Their relaxation is ensured 

' The lowest possible energy state for any photospheric magnetic field 
distribution which develops during the emergence process would locally (i.e. 
above the emergence region) con'espond to a potential field that has to contain 
one single nullpoint, owing to th e presence of a closed PIL embedded in a 
region of nearly vertical field (e.g. |Antiochos|1998j see Introduction). 
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Fig. 3. — Selected magnetic field lines outlining the first reconnection which occurs between the emerging flux rope and the outer coronal arcade. The bottom 
panels show a zoom into the emergence region. Contours of Bz > (B^ < 0) at the bottom plane {z = 0} are shown in white (black). Polarity inversion lines 
(PILs) are drawn in yellow. The negative polarity of the large-scale potential field is located on the left-hand side of the main PIL. All field lines are calculated 
from fixed footpoints on the left-hand side of the main PIL in all panels. From left to right: at onset of (t = 22), during (t = 34), and after (t = 43) the 
reconnection. Blue field lines the negative polarity of the emerging flux rope, green field lines show the core of the prominence flux rope. Red field lines outline 
the outer coronal arcade at early times and the newly formed small loop system later on. Pink field lines show the inner coronal arcade. ArTowheads mark the 
direction of the magnetic field. The black circles located above the main PIL mark field line dips, which are assumed to carry prominence material. (An mpeg 
animation of this figure is available in the online journal.) 



through the launch of a torsional Alfven wave which travels 
from low altitudes all along the arcade (see the evolution of 
the blue field lines in Figs.[3]and|4ji. 

Such reconnection between twisted and untwisted coro- 
nal fields has been suggested by several authors as a driving 
mechanism for jets (see, e.g., Schmieder et al.|1995 Shibata 



et al.||1997 l. If the reconnection is sufficiently impulsive, it 



can launch a shear (in 2.5D) or torsional (in 3D) Alfven wave, 
which can accelerat e the plasma up ward, as shown in numer- 
ical s i mulati ons by Yokoyama & Shibata (1999) and Pariat 
|et al.| ( |2009| l, respectively. ITie impulsive nature and large 
wavenumber of the wave in our simulation is a priori not ex- 
pected, since the transition from nearly unsheared to highly 
sheared emerging flux rope fields lines involved in the recon- 
nection is continuous. Also, our code does not incorporate 
any time-varying resistivity. Therefore, the impulsive launch 
of the wave must result from some perturbation of the system 
which yields a strong increase of the reconnection rate. In- 
deed, we find a strong expansion of the not yet reconnected 
central part of the emerging flux rope at i « 40 (i.e. between 
the stages shown in the middle and right panels of Fig. [3] see 
the corresponding online animation). At this time, a flux rope 



twist of 1.5 turns has entered the corona, indicating that 
this sudden increase in expa nsion might be relate d to the on- 
set of a kink instability (as in F an & Gibson|2003 2007). The 
reconnection-driven torsional Alfven wave in our simulation 
suggests an explanation for the jet-like brightening traveling 
along the cavity loops observed by XRT in our event, as well 
as for the observed transverse oscillation of the upper parts of 
the cavity loops (see Sect. 2.2 1. 

The transverse deformation of the blue field lines in our 
simulation during the passage of the wave (Fig.|4|l is, however, 
obviously much larger than what is observed in our eve nt and 
typically in coronal jet-like events ( Cirtain et al.|2007 i. This 
might simply be due to the fact that, although we have chosen 
the emerging flux rope to be as small as possible as compared 
to the prominence flux rope (see Sect.|3]l, the difference in size 
between the emerging flux and the prominence-cavity system 
might still be significantly larger. Also, to some extent the 
unrealistic reconnection time scales in our simulation might 
play a role. They are mostly constrained by the intrinsic dif- 
fusivity of the numerical scheme and by the prescribed mag- 
netic field smoothing, and do neither correspond to fast re- 
connection nor to the seminal Sweet-Parker regime. Still, the 
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Fig. 4. — Snapshots of the configuration outlining the second reconnection, which occurs between the newly formed small loop system and the inner coronal 
arcade. Field lines, contours, and black circles are as in Fig.ls] The bottom panels show a zoom into the emergence region. From left to right: at onset of (t = 46) 
and during (t = 49, 57) the reconnection. Pink field lines snow the inner arcade initially, and the second newly formed loop system later on. The propagation 
of an Alfven wave along the blue field lines and the slipping motion of the footpoints of the red field lines on the right-hand side of the prominence flux rope are 
visible. (An mpeg animation of this figure is available in the online journal.) 



qualitative agreement which we find here with other coronal 
jets simulations performed with different codes (see Sect. [TJ 
suggests that we can believe in the overall mechanism for jet 
acceleration which our simulation finds. 

Note that the nullpoint-related elbow in the lower part of the 
blue field lines apparently disappears during the evolution in 
our simulation (as it does in the observation too; see Fig. |2]i. 
In the simulation, this is merely due to projection effects and 
the motion of the nuUpoint. From a different viewing angle, 
an elbow at low heights remains visible. 

Up to this point, the evolution is as expected from the clas- 
sical model for coronal jets and previous 2.5D and 3D simu- 
lations of this model (see Sect.[T]and Fig.[5^,b). It results in 
the formation of a 3D nullpoint topology, but actually half of 
its fan surface still consists of not reconnected emerged field 
lines, so that only the other half is expected to brighten in soft 
X-rays. In other words, only a /laZ/ anemone has been formed 
at this stage in the simulation. 

4.2. Second step: formation of the second half of the 
anemone 

As the reconnection-driven transfer of twist/shear from the 
emerging flux rope into the ambient coronal field progresses, 
the footpoints of the newly formed red field lines located 
within the positive polarity of the emerging flux rope are con- 



tinuously displaced in the negative y direction (see Fig.[3]and 
the online animation), toward inner arcade field lines which 
also overlie the prominence flux rope and are yet unaffected 
by the reconnection (shown in pink in Fig. 3}. Meanwhile, 
following the magnetic field reorientation at nigh altitude, the 
nullpoint (and hence the reconnection site) undergoes a coun- 
terclockwise horizontal rotation, from the leftmost edge of the 
emerging rope toward its center, thereby approaching the pink 
arcade field lines (see Fig. |5]). The displacement of the foot- 
points of the small red loop system corresponds to the appar- 
ent expansion of LI toward L2 in the observation, which starts 
at about 18:45 UT (see Sect.[Z2l). 

Eventually a second reconnection starts, now between the 
previously reconnected sheared red field lines and the pink 
arcade field lines (Fig. Hli. It must, and indeed does, take 
place at the nullpoint, which has gradually moved toward the 
pink arcade during the first reconnection episode (see above). 
This second reconnection leads to the formation and growth 
of a second small loop system (cusp-shaped pink field lines 
in Fig. 4|i, which corresponds to the growth of the observed 
bright loop system L2 from about 18:50 UT on (Fig.|2]i. The 
reconnected red and pink loop systems both have footpoints 
within the positive polarity of the emerging flux, and they col- 
lectively display an anemone-like shape which is significantly 
wider than this polarity (Fig. [5};; compare also with the col- 
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Fig. 5. — (a-c): Side view on tlie flux emergence region in tlie simulation at t 
Selected magnetic field lines are shown in order to outline the topological evolution, (d): Top view on the center of the anemone-shaped loop pattern i 
(c), outlining the twist of the field lines, (e) and (f): Top and side view on the emergence region at t = 51, showing field lines starting from a circle around the 
spine (black line) with twice the radius as used in (c) and (d). These field lines approximately outline the fan surface. The respective positions of the nuUpoint 
are indicated by circles. 



17, 34, and 51. The contours at the bottom plane are as in Figs.[3]and|4] 

I snown in 



lective shape of LI and L2 in Fig.|2p,e). 

In the course of the second reconnection, the sheared red 
loop system "feeds" the newly formed small pink loop sys- 
tem with magnetic flux, hence the former shrinks while the 
latter grows (Fig. HI). Although, at first glance, the bottom 
right panel gives the impression that the red loop system has 
almost disappeared, the elbow visible in the rightmost red ar- 
cade field line shows that there still exist relatively extended, 
not reconnected red loop field lines at this late phase of the 
second reconnection. These field lines are rooted between the 
elbow-shaped field line and the rightmost, very flat red field 
line. The final relative size of the red and pink loop systems 
will depend on how long reconnection persists, which, in turn, 
will depend on a number of factors, as for example on the ratio 
of the respective magnetic fluxes present in the reconnecting 
flux systems and on the lowest state of magnetic energy the 
system can reach by means of its relaxation. A shrinking of 
LI during the growth of L2 is also indicated by the obser- 
vations (compare Fig. |2|l and e). However, this "shrinking" 
might be merely due to a relatively fast cooling of LL Study- 
ing the issues of lifetime and visibility in soft X-ray of the 
reconnected loops in our simulation would require the inclu- 
sion of proper thermodynamics, which is beyond the scope of 
this paper. 



4.3. Geometrical properties of the reconnected field 

The geometrical properties of these new loop systems in 
the simulation are different from what has been found in the 
simulations mentioned in Sect.[T| They arise from a two-step 
and fully 3D transfer of sheared flux, first from the sheared 
core of the emerging flux rope into the new small loop system 
formed at its side (as the emerging blue field lines reconnect 
with ambient [red] arcade field lines that are anchored on the 
left-hand side of the emerging rope), and second from the very 
same loop system into a second generation of small loops, 
which form at the other side of the emerging flux rope (as 
the sheared small red field lines reconnect with other ambient 
[pink] arcade field lines). 

As a result, a full anemone forms around the parasitic po- 
larity of the emerging bipole. A large fraction of the nullpoint 
associated fan now consists of once or twice reconnected field 
lines, and not to a large extent of emerged field lines as in most 
of the models described in Sect.[T] Since both sides of the fan 
have been formed through magnetic reconnection, they can a 
priori be both visible as hot loops in soft X-rays. The center 
of the anemone structure contains significant twist once it has 
fully formed, although the structure appears to be potential 
when viewed from some distance (Fig. |5t). The twist is con- 
centrated around the inner spine of the nullpoint (see Fig.jSH), 
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but some of it is also present along the fan. This twist is the 
remnant of magnetic shear that has not been ejected in the 
form of a torsional Alfven wave along the large-scale recon- 
nected arcades that overlie the prominence flux rope. 

It is worth noting that, even though a fully 3D anemone has 
been formed, it still locally possesses a quasi, but not exact, 
translational invariance along the y axis, i.e. along the axis 
of the emerging flux rope, around the nullpoint. This can be 
seen in Fig. l5b,d: almost all of the red and pink field lines 
fan out roughly along the x direction, i.e. perpendicular to 
the rope axis, even though they have been integrated from 
equidistant footpoint positions along a small circle centered 
at the inner spine of the nullpoint, within the positive polarity 
of the emerging flux, where is roughly constant (as shown 
in Fig. ISH). This strong departure from axisymmetry around 
the nullpoin t, which does not exist in afl of the 3D nullpoint 
models (e.g. Antiochos'1998VPariat et al."2009\ can here be 
attributed to the elongation of the emerging flux, and possibly 
to the different inclinations of the ambient field lines. The few 
red and pink field lines in Fig.|5|J which fan out from the null- 
point region along the positive y direction, clearly show that 
the anemone does not contain a true null line, but one sin- 
gle nullpoint. The quasi-translational invariance is then due 
to fan-related eigenvalues of the single nullpoint, which have 



very different amplitudes ( Lau & Finn|1990^ . This property, 
which has already been identified in the simulation of Masson 
et al.| (2009 ) with a different MHD code, may be present also 
in other simulations, such as those w hich let emerge a very 
long flux rope (e.g. More no-Insertis et al.|2008| . 

Fig. and f show that, despite the strong departure from 
axisymmetry, the fan surface also extends in y-direction, i.e., 
along the axis of the emerging flux rope. The field lines shown 
in the two panels are located very close to the fan surface and 
outline a "heart-shaped" form of the fan. We expect the fan 
to develop a more uniform radial distribution as the system 
relaxes to a force-free field, or if less el ongated parasitic po- 
larities are considered ( |Antiochos| 1 998| l. 

4.4. Response of the prominence cavity to the anemone/jet 
formation 

While the two-step reconnection described in the previous 
subsections persists beside the prominence flux rope, the tor- 
sional Alfven wave triggered during the first reconnection 
phase travels upward (see the evolution of the blue field lines 
in the top panels of Fig.[4]|. Even though the upper part of the 
prominence flux rope significantly bends to the side while the 
wave passes by, the simulated "prominence material" located 
in the flux rope dips (compute d as in [A ulanier & Demoulin 
T998| lyan Ballegooijen|[2004l |Gibson & Fan 2006, .Dud& | 
et al.|2008| and indicated by the black circles in Figs.|3[and|4| 
does not show a significant motion. This difference is due to 
the relatively small twist of the rope combined with its non- 
negligible curvature, which result in the occurrence of mag- 
netic dips only far below the rope axis, at altitudes low enough 
not to be significantly affected by the wave (see also the dis- 
cussion onthe strength of the perturbation caused by the wave 
in Sect.pTT]). 

A perturbation of the cavity is observed by XRT on the 
right-hand side of the prominence in our event. Between 
19:01 and 19:19 UT, after the formation of LI, while the jet- 
like brightening is still propagating and L2 is still growing, 
a continuous motion ("slipping") of the lowermost parts of 
several loops, from the outer edge of the cavity toward the 
prominence, becomes visible (see Fig.|2^, Fig.jSJ", and the cor- 




FlG. 6. — Top and oblique view on two field lines within the emerging 
flux rope alt = 31, outlining the change of field line inclination during the 
emergence. See the text for details. 



responding XRT movie). Estimating their slipping velocities 
by following them individually is not straightforward, as all 
these loops are not very much contrasted with respect to the 
background corona. The only sure number we could derive 
is a minimum drift velocity of 35 kms~^ for the ensemble of 
these loops, since they all have moved by 51" along the solar 
limb during a time interval of 1 8 minutes. Still, we cannot rule 
out that the drift velocity of individual loops could be much 
larger. 

The very same phenomenon takes place in the simulation. 
The online animation of the simulation shows a slipping mo- 
tion of the footpoints of the red large-scale arcade field lines 
on the right-hand side of the prominence flux rope (see also 
the top panels in Fig.[4|l. This slipping occurs along the foot- 
points of arcade field lines (of which some are shown by the 
pink field lines in the left panels of Fig. Hli. It starts right af- 
ter the second reconnection sets in and begins to exchange 
the connectivities between the small (and already once re- 
connected) sheared red field lines and large-scale arcade field 
lines. This slipping of field lines, after they have reconnected 
at the nullpoint, is not a numerical artifact, caused by insuffi- 
cient accuracy of the field line integration or by a too strong 
smoothing of the magnetic field during the MHD simulation: 
the slipping occurs over many grid points in the large-scale 
positive polarity beside the prominence flux rope, and it does 
not occur outside of flux regions which reconnect at the null- 
point. A very mild slipping is noticeable also during the first 
reconnection, at the footpoints of the reconnected large-scale 
blue lines along which the wave travels (see the online anima- 
tion). 

The explanation of t his sUp ping is very p robably the same 
as first put forward by Masson et al. ( 2009) in their nullpoint 
associated confined flare model: the asymmetry of the fan- 
spine configuration, manifesting as a local quasi-translational 
invariance of the magnetic field around a nullpoint, which has 
one fan-r elated eigenvalue of very small, but finite value (see 
Sect. |4.2| and Figs.|5]:,d), results in the appearance of a narrow 
halo of finite width around the nullpoint separatrices, in which 
field lines have strong squashing de grees, i.e. which consti - 
tutes quasi-separatrix layers (QSLs; DemouUn et al. 19961. 
It has been recently shown that current sheet formation and 
diffusion at QSLs result in slipping/slip-running magnetic 
reconnection, manifesting as apparent field line motions at 
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sub/super Alfvenic speeds ( Aulanier et al. 2005 2006 1. In the 
case of a 3D nullpoint topology, the presence of a QSL halo 
surrounding both the fan surface and the spine field line re- 
sults in a complex pattern for magnetic reconnection, in which 
a given field line slips and slip-runs, both before and after its 
connectivity jumps at the nullpoint. To the best of ou r knowl- 
edge, the pr esent simulation is the second one after Masson 
et al. ( 2009 ) for which this sequential nullpoint and slipping 
reconnection behavior is reported. Based on these simula- 
tions, the slipping loops observed by XRT in our event can be 
regarded as a new direct evidence for slipping reconnection in 
the corona, as previously observed by |Aulanier et al.| ( |2007| l. 

5. DISCUSSION 

We presented a /3 = 3D MHD simulation of the interac- 
tion of a small emerging bipole with a large-scale arcade-like 
coronal field overlying a weakly twisted coronal flux rope. 
The simulation was developed in order to understand Hin- 
odelXKJ observations of a small limb event which took place 
at the edge of a quiet Sun prominence cavity. The event 
showed a number of puzzling features, which could satisfac- 
torily be explained by the magnetic field evolution in the sim- 
ulation (see Sect.[4]for details). 

The main results of the simulation can be summarized as 
follows: 

1. The emergence of a twisted flux rope into one large- 
scale polarity of an arcade-like coronal field yields the 
formation of a fully 3D single nullpoint topology in the 
corona, consisting of a fan-spine configuration in which 
the fan surface significantly extends over the parasitic 
polarity of the emerged flux. 

2. The configuration forms in a two-step reconnection pro- 
cess at the nullpoint, which yields the successive for- 
mation of two small loop systems below the fan sur- 
face, on opposing sides of the parasitic polarity. The 
first loop system hereby "feeds" the second one with 
magnetic flux. Since the two loop systems (and field 
lines surrounding them) have common footpoints in the 
parasitic polarity, they collectively display an anemone 
shape. 

3. The reconnection facilitates the transfer of twist/shear 
from the emerging flux rope into the coronal arcades by 
means of a torsional Alfven wave which travels along 
the arcades. A fraction of the twist remains below the 
fan surface, where therefore the resulting magnetic field 
is nonpotential. 

4. The wave is launched by reconnection, as the expansion 
velocities of the emerging fields in the corona suddenly 
increase, possibly due to the onset of a ideal MHD 
kink instability in the not yet reconnected part of the 
emerging flux rope, after sufficient twist has entered the 
corona. 

5. The 3D nullpoint reconnection is accompanied by slip- 
ping reconnection of arcade field lines on the oppo- 
site side of the pre-existing coronal flux rope. This 
can probably be explained by the presence of quasi- 
separatrix layers around the nullpoint. 

Our simulation combines dynamic and topological ele- 
ments of flux emergence into unipolar coronal fields and coro- 



nal jet formation which have previously not been found alto- 
gether in a single simulation. The two-step reconnection pro- 
cess found in our simulation provides a new model for the for- 
mation of 3D nullpoint topologies with extended fan surfaces 
in the corona. By means of reconnection between twisted and 
untwisted fi elds, the model ca n also acc ount for coronal jet 
activity (e.g . Schmieder et al. 1995 ; Shibata et al. 1997[ Pariat 
et al.||2009 l. Finally, since the fields below the fan surface 
are all fornied by reconnection, and can therefore a priori as- 
sumed to brighten in soft X-rays, our simulation also suggests 
a mechanism for the formation of (full) anemone active re- 
gions. The fact that some of the emerging twist remains be- 
low the fan surface in the final configuration supports models 
of jet formation which assume a twisted n ull point configura- 
tion prior to the jet (e.g. 'P ariat et al.|2009j l. 

Our model can, however, not account for the long-lived evo- 
lution of anemone active regions. Such regions are observed 
to consist of bright loops on both sides of a paras itic polar- 
ity over time-periods of days (e.g. Asai et al.|2008 1, whereas 
the reconnection in our simulation is expected to produce 
rather short-lived brightenings, which would occur succes- 
sively rather than simultaneously (as in the event described 
in this paper). It seems likely that the long-lived dynamics of 
anemone regions is due to continuous dynamic perturbations 
of the configuration after its formation by, for instance, ongo- 
ing flux emergence and photospheric flows. Our model might 
serve as a starting point for a numerical study of anemone re- 
gion dynamics. 

Let us now briefly discuss the robustness of our results with 
respect to the variation of various model parameters. A de- 
tailed investigation of this question will require an extensive 
parametric study, but some aspects may be discussed here. 

First of all, the model relies on fully 3D effects. The most 
important ingredient is the twist brought up by the emerg- 
ing flux rope. For the two-step reconnection as described in 
this paper to work, the emerging field lines must significantly 
change their inclination as they rise into the corona, such that 
the first reconnected loop system can evolve toward yet un- 
reconnected remote ambient fields during the reconnection. 
If the twist in the emerging flux rope (more precisely: the 
change of field line inclination along the cross-section of the 
rope) is chosen too small, the first reconnection will occur 
(as in the flux emergence simulations mentioned in Sect. [T]), 
but the second one will not set in. Our adopted boundary- 
driven "kinematic" flux emergence technique models the rigid 
emergence of a twisted flux rope into the corona, and does 
therefore not capt ure the d ynamics of real flux emergence (for 
a discussion, see |Fan & Gibson 2003 ). In particular, it as- 



sures that the inclmation of the emerging field lines is contin- 
uously changing. It would be, of course, desirable to test our 
m odel in a "dynamic" flux emer gence simulation, as the one 
by [More no-Insertis et al.j ( |2008[ ). Here we have chosen the 
"kinematic" approach since, at present, "dynamic" flux emer- 
gence simulations do not allow to consider complex coronal 
magnetic field configurations as the one used here, which was 
chosen to model the observed prominence-cavity system. 

Apart from the twist profile of the emerging flux rope, there 
are several other model parameters which had to be chosen 
carefully in order to match the observations (but not neces- 
sarily to produce a 3D nullpoint configuration, which is more 
general). First, of course, the size of the emerging flux rope 
and its distance from the coronal flux rope had to be adjusted 
as suggested by the observations. Second, the choice of the 
magnetic field strength within the emerging rope (and there- 
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fore its magnetic flux ) is important. The ratio of the field 
strength (and flux) between the emerging rope and the ambi- 
ent corona will influence, for instance, the amplitude of the 
torsional Alfven wave and the final size of the reconnected 
small loops, i.e. the size of the anemone. In the simulation 
presented here, the field strength within the coronal part of 
the emerging flux rope was typically about five times larger 
than in the neighboring arcade fields throughout most of the 
evolution, so that the flux within the (small) rope was not neg- 
ligible compared to the flux contained in the large-scale ar- 
cades. Third, the inclination of the initial coronal field might 
play a role, although we expect a similar evolution if the flux 
rope would be emerged into a purely vertical coronal field, 
mimicking "open" field lines in coronal holes. Finally, the 
magnetic orientation of the emerging rope with respect to the 
pre-existing coronal configuration plays an important role. If 
this quantity would be reversed in the simulation, and if the 
evolution would be seen from the same viewing angle as in 
Figs. |3] and H] we expect the second reconnected (pink) loop 
system to form on the left-hand side of the first (red) one, 
and the blue reconnected field lines above the nullpoint to de- 
velop an elbow which bends away from the prominence flux 
rope rather than toward it. It would therefore be interesting to 
study how the system behaves if intermediate orientations of 
the emerging rope are chosen. 

Our results underline the importance of a precise examina- 
tion of the magnetic topology (and of its formation) for the 
understanding of many dynamic events in the solar corona. 



Antiochos, S. K. 1998, ApJ, 502, L181 

Archontis, V., Galsgaard, K., Moreno-Insertis, R, & Hood, A. W. 2006, ApJ, 
645, L161 

Asai, A., Shibata, K., Hara, H., & Nitta, N. V. 2008, ApJ, 673, 1 188 
Aulanier, G., DeLuca, E. E., Antiochos, S. K., McMullen, R. A., & Golub, 

L. 2000, ApJ, 540, 1126 
Aulanier, G., & Demoulin, P. 1998, A&A, 329, 1125 
Aulanier, G., Demoulin, P., & Grappin, R. 2005, A&A, 430, 1067 
Aulanier, G., Golub, L., DeLuca, E. E., Cirtain, J. W., Kano, R., Lundquist, 

L. L., Narukage, N., Sakao, T., & Weber, M. A. 2007, Science, 318, 1588 
Aulanier, G., Paiiat, E., Demoulin, P, & Devore, C. R. 2006, Sol. Phys., 

238, 347 

Bateman, G. 1978, MHD instabilities (Cambridge, Mass., MIT Press) 
Cirtain, J. W., Golub, L., Lundquist, L., van Ballegooijen. A., Savcheva, A., 
Shimojo, M., DeLuca, E., Tsuneta, S., Sakao, T., Reeves, K., Weber, M., 
Kano, R., Narukage, N., & Shibasaki, K. 2007, Science, 318, 1580 
Craig, I. J. D., & McClymont, A. N. 1991, ApJ, 371, L41 
Demoulin, R, Henoux, J. C, Priest, E. R., & Mandrini, C. H. 1996, A&A, 
308, 643 

Demoulin, P., van Driel-Gesztelyi, L., Schmieder, B., Hemoux, J. C, 

Csepura, G., & Hagyard, M. J. 1993, A&A, 271, 292 
Dudi'k, J., Aulanier, G., Schmieder, B., Bommier, V., & Roudier, T. 2008, 

Sol. Phys., 248, 29 
Fan, Y., & Gibson, S. E. 2003, ApJ, 589, L105 
— . 2007, ApJ, 668, 1232 

Feynman, J., & Martin, S. F 1995, J. Geophys. Res., 100, 3355 
Filippov, B. 1999, Sol. Phys., 185, 297 

Filippov, B., Golub, L., & Koutchmy, S. 2009, Sol. Phys., 254, 259 

Fletcher, L., Metcalf, T. R., Alexander, D., Brown, D. S., & Ryder, L. A. 
2001, ApJ, 554, 451 

Gibson, S. E., & Fan, Y. 2006, Journal of Geophysical Research (Space 
Physics), 111, 12103 

Golub, L., Deluca, E., Austin, G., Bookbinder, J., Caldwell, D., Cheimets, 
P, Cirtain, J., Cosmo, M., Reid, P, Sette, A., Weber, M., Sakao, T., Kano, 
R., Shibasaki, K., Hara, H., Tsuneta, S., Kumagai, K., Tamura, T., 
Shimojo, M., McCracken, J., Carpenter, J., Haight, H., Siler, R., Wright, 
E., Tucker, J., Rudedge, H., Barbera, M., Peres, G., & Varisco, S. 2007, 
Sol. Phys., 243, 63 

Gombosi, T. 1., Powell, K. G., & de Zeeuw, D. L. 1994, J. Geophys. Res., 
99, 21525 



Without a detailed consideration of the topology, it would 
have been very difficult to understand the complex sequence 
of dynamic features occurring in the simulation and in the ob- 
served event. 



We thank the referee for detailed comments and suggestions 
which helped very much to improve the quality of this paper. 
Financial support by the European Commission through the 
SOLAIRE Network (MTRN-CT-2006-035484) and through 
the FP7 SOTERIA project (Grant Agreement n° 218816) are 
gratefully acknowledged. Hinode is a Japanese mission de- 
veloped and launched by ISAS/JAXA, with NAOJ as domes- 
tic partner and NASA and STFC (UK) as international part- 
ners. It is operated by these agencies in co-operation with 
ESA and the NSC (Norway). LG and KKR are supported by 
NASA under contract NNM07AB07C to SAO. We thank all 
the people who collaborate actively within the JOP 178 ob- 
servations campaign, particularly T. Berger and the Hinode 
team, G. Molodij and the Meudon Solar Tower observers, and 
P. Mein for reducing the MSDP data, and T. Roudier and N. 
Labrosse for the coordination and the pointing of the observa- 
tions. TT thanks Y. Fan for several discussions on "kinematic 
flux emergence". We also thank P. Demoulin, B. Kliem, and 
E. Pariat for very helpful comments on this work. Our MHD 
calculations were done on the dual-core quadri-Opteron com- 
puters of the Service Informatique de I'Observatoire de Paris 
(SIO). 



Heinzel, P., Schmieder, B., Famfk, R, Schwartz, P., Labrosse, N., Kotrc, P., 
Anzer, U., Molodij, G., Berlicki, A., DeLuca, E. E., Golub, L., Watanabe, 
T, & Berger, T. 2008, ApJ, 686, 1383 

Heyvaerts, J., Priest, E. R., & Rust, D. M. 1977, ApJ, 216, 123 

Isenberg, P A., & Forbes, T. G. 2007, ApJ, 670, 1453 

Karpen, J. T., Antiochos, S. K., Devore, C. R., & Golub, L. 1998, ApJ, 495, 
491 

Keppens, R., Nool, M., Toth, G., & Goedbloed, J. R 2003, Computer 

Physics Communications, 153, 317 
Kliem, B., Karlicky, M., & Benz, A. O. 2000, A&A, 360, 715 
Kliem, B., & Torok, T. 2006, Phys. Rev. Lett., 96, 255002 
Kosugi, T, Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., Tachikawa, S., 
Hashimoto, T, Minesugi, K., Ohnishi, A., Yamada, T., Tsuneta, S., Hai'a, 
H., Ichimoto, K., Suematsu, Y, Shimojo, M., Watanabe. T., Shimada, S., 
Davis, J. M., Hill, L. D., Owens, J. K., Title, A. M., Culhane, J. L., Harra, 
L. K., Doschek, G. A., & Golub, L. 2007, Sol. Phys., 243, 3 
Lau, Y.-T., & Finn, J. M. 1990, ApJ, 350, 672 
Low, B. C, & Hundhausen, J. R. 1995, ApJ, 443, 818 
Masson, S., Pariat, E., Aulanier, G., & Schrijver, C. J. 2009, ApJ, 700, 559 
Moreno-Insertis, R, Galsgaard, K., & Ugarte-Urra, I. 2008, ApJ, 673, L211 
Murray, M. J., van Driel-Gesztelyi, L., & Baker, D. 2009, A&A, 494, 329 
Nishizuka, N., Shimizu, M., Nakamura, T., Otsuji, K., Okamoto, T. J., 

Katsukawa, Y, & Shibata, K. 2008, ApJ, 683, L83 
Nistico, G., Bothmer, V., Patsourakos, S., & Zimbardo, G. 2009, ArXiv 
e-prints 

Pariat, E., Antiochos, S. K., & DeVore, C. R. 2009, ApJ, 691, 61 
Patsourakos, S., Pariat, E., Vourlidas, A., Antiochos, S. K., & Wuelser, J. P. 

2008, ApJ, 680, L73 
Pontin, D. I., Bhattacharjee, A., & Galsgaard, K. 2007, Physics of Plasmas, 

14, 052106 

Priest, E. R., & Forbes, T. G. 2002, A&A Rev., 10, 313 

Sato, T., & Hayashi, T. 1979, Physics of Fluids, 22, 1 189 

Schmieder, B., Shibata, K., van Driel-Gesztelyi, L., & Freeland, S. 1995, 

Sol. Phys., 156, 245 
Schrijver, C. J., Elmore, C, Kliem, B., Torok. T., & Title, A. M. 2008, ApJ, 

674, 586 

Shibata, K., Ishido, Y, Acton, L. W., Strong, K. T., Hirayama, T, Uchida, 
Y, McAllister, A. H., Matsumoto, R., Tsuneta, S., Shimizu, T., Hara, H., 
Sakurai, T., Ichimoto, K., Nishino, Y, & Ogawara, Y. 1992a, PASJ, 44, 
L173 



12 



Torok et al. 



Shibata, K., Nitta, N., Strong, K. T., Matsumoto, R., Yokoyama, T., 
Hirayama, T., Hudson, H., & Ogawara, Y. 1994, ApJ, 431, L51 

Shibata, K., Nozawa, S., & Matsumoto, R. 1992b, PASJ, 44, 265 

Shibata, K., Shimojo, M., Yokoyama, T., & Ohyama, M. 1997, in 
Astronomical Society of the Pacific Conference Series, Vol. Ill, 
Magnetic Reconnection in the Solar Atmosphere, ed. R. D. Bentley & 
J. T. Mariska, 29-38 

Titov, V. S., & D^mouUn, R 1999, A&A, 351, 707 

Torok, T., & Kliem, B. 2003, A&A, 406, 1043 

— . 2005, ApJ, 630, L97 



— . 2007, Astron. Nachr., 328, 743 
Torok, T, Kliem, B., & Titov, V. S. 2004, A&A, 413, L27 
Ugarte-Urra, I., Warren, H. R, & Winebarger, A. R. 2007, ApJ, 662, 1293 
van Ballegooijen, A. A. 2004, ApJ, 612, 519 
Yokoyama, T, & Shibata, K. 1996, PASJ, 48, 353 
Yokoyama, T, & Shibata, K. 1999, in American Institute of Physics 
Conference Series, Vol. 471, American Institute of Physics Conference 
Series, ed. S. T. Suess, G. A. Gary, & S. F. Nemey, 61-66 



